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In many braneworld models, gravity is largely modified at the electro-weak scale ~ ITeV. In 
such models, primordial black holes (PBHs) with a lunar mass M ~ 1O" 7 M might have been 
produced when the temperature of the universe was at ~lTeV. If a significant fraction of the dark 
halo of our galaxy consists of these lunar mass PBHs, a huge number of BH binaries will exist in our 
neighborhood. Third generation detectors such as EURO can detect gravitational waves from these 
binaries, and can also determine their chirp mass. With a new detector designed to be sensitive at 
high frequency bands > 1kHz, the existence of extradimensions could be confirmed. 



Introduction) In recent years, there has been a great 
interest in braneworld scenarios in which the ordinary 
matter fields are confined in a four dimensional object 
( "brane" ) , while the gravitational field propagates in the 
higher dimensional spacetime ( "bulk" ) . One motivation 
for the brane world scenario is solving the hierarchy prob- 
lem between the Planck scale and the electro- weak scale 
(~ ITeV). The possibility of relating the Planck scale to 
the TeV scale by considering large sized extradimensions 
was pointed out [1]. In the model proposed by Ran- 
dall and Sundrum (RSII) [2] one of the extradimensions 
extends infinitely, but compactification of the extradi- 
mension is effectively achieved by the warped geometry. 
In this model the four dimensional Planck mass M p is 
related to the brane tension a and the bulk curvature 
length I as M 2 = 47rtr£ 2 /3. Hence, keeping the scale 
of the brane tension at the TeV scale, the Planck scale 
can be derived by setting I w 0.1mm, although this fact 
does not mean that the RSII model solves the hierarchy 
problem. 

In such braneworld scenarios, the evolution of the uni- 
verse at a temperature above ~ ITeV can be dramati- 
cally altered. For instance, the brane universe may have 
had an era of violent "mesoscopic" activity at scale as 
large as I ~ 1mm at a temperature T > ITeV [3]. Such 
a violent activity may generate large-amplitude fluctu- 
ations at the horizon scale H~ l . Then a large num- 
ber of primordial black holes (PBHs) with lunar masses 
M B h ~ H~ l Ml ~ 10 -7 M© could be produced. To date, 
PBH physics in the context of braneworld scenarios has 
been studied [4,5], but none of the literature has shed 
any light on the astrophysical implication of considering 
PBHs as cold dark matter. 

Currently, there is no stringent observational con- 
straint on SUb-Lunar mass Compact Objects (SULCOs) 
with mass 1O~ 12 M < M < 10~ 7 Af o as the dark matter 
candidates, which can be either molecular clouds, small 
planets, or PBHs. Microlensing tests such as the EROS 
and MACHO collaborations ruled out the possibility that 
compact objects with mass 10~ 7 M S < M < 10~ 3 M Q 
make up the halo dark matter in the Milky Way [6]. 



On the other hand, femtolensing tests of Gamma-ray 
bursts (GRBs) ruled out the mass range 10 _16 Mq < 
M < 10 _13 Mq in the universe assuming that GRBs are 
at cosmological distance so that the angular size of the 
GRB source is sufficiently small [7], SULCOs also in- 
duce picolensing of GRBs, but the observational limit is 
very weak [7]. Dynamical constraints on the number of 
SULCOs in the dark halo are also less stringent [8] . 

In this letter, we consider the PBH SULCOs that have 
been produced by violent activity at temperatures ~TeV 
as the dominant constituent of the cold dark matter. 
Then there must be a huge number of PBHs (> 10 19 ) 
in the dark halo of the Milky Way. As is discussed in 
the PBH MACHO scenario [9], it is natural to expect 
that most of the PBHs form binaries and some of them 
coalesce emitting gravitational waves within the present 
age of the universe. From the frequency evolution in such 
gravitational waves, we may be able to probe the devi- 
ation from Newtonian gravity typical of scenarios with 
infinite extradimensions. 

Formation and evolution of PBH SULCO binaries) We 
consider the formation of PBHs in the context of grav- 
ity theories which have a critical energy scale at the 
TeV scale. If PBHs are formed at a temperature T c ~ 
TeV, they are of the lunar mass, Mbh ~ H~ 1 Mp ~ 

0.1mm Ml {T c /TeVy 2 ~ 10- 7 M Q (T c /TeV) -2 . Here we 
have assumed that the number of effective degrees of free- 
dom at the PBH formation epoch is about 100. Although 
it has been argued that PBHs could undergo substantial 
growth by accreting material from the cosmological back- 
ground in the high-energy regime (p > a) [4,5], such an 
effect is not relevant in our discussion because here PBHs 
are formed right after the end of the high-energy regime. 

In order to estimate the event rate of coalescence of 
PBH SULCO binaries, we follow the arguments in [9]. 
For simplicity, we assume that the PBH SULCOs dom- 
inate the dark matter, i.e., Q m — £Ipbh, and the PBH 
mass spectrum is monochromatic. 

At matter-radiation equality, the comoving mean sep- 
aration is x = 2 x lO 16 (M sff /M ) 1 / 3 (fi m /i 2 )- 4 / 3 cm, 
where h is the Hubble parameter in units of 100km 



s _1 Mpc _1 . Henceforth, we adopt £l m h 2 — 0.15 as sug- 
gested by recent observations. A binary with its sepa- 
ration x decouples from the cosmic expansion when the 
local energy density becomes a few times larger than the 
energy density of the background radiation. The tidal 
force from the nearest PBH at a comoving distance y to 
the center of mass of the binary adds angular momentum 
to keep the holes from a head-on collision. For simplicity, 
we assume that the distribution of PBHs are spatially ho- 
mogeneous and isotropic. Then the distribution of x and 
y takes the form p(x,y) oc x 2 y 2 exp(— y 3 /x 3 ), which is 
further approximated by a uniform distribution oc x 2 y 2 
in the range x < y < x. 

The coalescence time t due to emission of gravitational 
waves can be written in terms of the semi-major axis a 
and the eccentricity e as t = to(a / do) 4 (1 — e 2 ) 7 ' 2 , where 
t Q = 10 10 yr and a = 2.4 x 10 11 (M BH /M Q ) 3 / 4 cm is the 
initial separation of a circular binary which coalesces in 
to. Thus, writing x and y in terms of t, e, and x, one 
obtains the distribution of the coalescing time: 
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where i max = £o(^/ a o) 4 - Let us denote the remaining 
time before coalescence by T. For small T, we can ne- 
glect the eccentricity e because the radiation reaction acts 
to reduce it. Then, T has a one-to-one correspondence 
with the scmimajor axis, and hence with the frequency of 
the emitted gravitational wave f g . The explicit relation 
becomes 
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Hence, provided that T <C t a , the probability distribu- 
tion for log f g is obtained as 
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Gravitational waves from PBH SULCO binaries) Now we 
discuss observability of gravitational waves from PBH 
SULCO binaries. We assume that the local dark halo 
density in the solar neighborhood is 0.0079M Q pc~ 3 [10]. 
Let the total number of PBHs within the distance D s be 
N(D S ). Then the distance D m i n such that we can expect 
at least one coalescence event within T at D < D m [ n is 
obtained by solving N(D m i n )f g P(f g ) = 1. Thus we have 
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After averaging over the orbital period and the orien- 
tations of the binary orbital plane, the amplitude of 



gravitational waves from a coalescing binary in a circu- 
lar orbit with distance D at frequency f g is h ~ 6 x 
10- 23 (D/2OMpc)- 1 (^B///M ) 5 /3(/ g /ioOHz) 2 / 3 . Sub- 
stituting D here with -D m ; n , we obtain the characteris- 
tic amplitude h c = h^/n, where n — f g AT is the num- 
ber of cycles during the observation time AT. Then the 
wave strength of gravitational waves from a monochro- 
matic source observed by using an interferometer h s = 
/i c /(5/ s ) 1 / 2 is evaluated as 
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where we have taken into account the antenna pattern 
of the detector sensitivity averaged over the direction to 
the binary by multiplying the factor l/y/E. The rms of 
the signal to noise ratio averaged over the source direc- 

~ 1/2 

tions and orientations is given by h s (f g )/S n (f g ), where 
Sn 2 (f) is the strain sensitivity in unit of Hz -1 / 2 . Note 
that the above estimate for h s is valid only for frequen- 
cies below f s = 3 x 10 3 (M BH /M Q )- 5 / 8 (T/5 yr)- 3 / 8 for 
which the source frequency does not change rapidly dur- 
ing the observation time AT <C T. 

The planned spectral noise density of the European 
Gravitational Wave Observatory (EURO) is [11], 

S n (f) = 1(T 50 x [(//245Hz)- 4 + (./7360Hz)- 2 

+ (A/770Hz)(l+ f 2 // 2 )]Hz-\ 

where fk is the knee-frequency. As seen from Fig.l, 
gravitational waves from PBH SULCO binaries are de- 
tectable by such an interferometer for integration time 
AT > 5yr. Furthermore, we can measure the chirp mass 
M c = (M 1 M 2 ) 3/5 /(M 1 + M 2 ) 1/5 of each binary, where 
Mi and M 2 denote the masses of respective PBHs. This 
is because the change in frequency during the observa- 
tion time AT, Av = (df g /dt)AT, becomes sufficiently 
large compared with the frequency resolution^ 1/AT 
at f g 3> 2Hz for AT = 5yr. Thus, the existence of 
SULCO binaries can be confirmed by observing gravi- 
tational waves. 

Probing extradimensions from gravitational waves) If 
PBH SULCO binaries are detected, we will have a chance 
to probe the existence of extradimensions. Here we fo- 
cus on the RSII model as a typical example in which 
the Kaluza-Klcin (KK) spectrum is continuous with no 
mass gap. In our present context, the brane tension 
is to be set to ^(TcV) 4 since it is expected to be re- 
lated to the energy scale of the Standard Model fields 
localized on the brane. Then, the observed four dimen- 
sional Planck mass is reproduced by setting the bulk 
curvature to the observationally allowed maximum value 
~ 0.1mm. The cosmic expansion law in this model is 
given by H 2 = &ir(p + 2p 2 /a)/3M 2 , where p is the en- 
ergy density of the universe. When p > a, the cosmic 
expansion law dramatically changes. 

Let us consider the non-relativistic gravitational po- 
tential V for a binary that consists of two point particles 



2 



with mass M 1 and M 2 on the brane. Let M — Mi + M 2 
be the total mass and /i = M1M2/M be the reduced 
mass. The separation between the two particles is de- 
noted by r. Continuous KK modes add a correction to 
the standard Newtonian term as [12], 
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provided that the separation r is much larger than the 
AdS radius I. Therefore, the correction owing to the con- 
tinuous KK modes is second order in e = GM/r, i.e., the 
second post-Newtonian (2PN) order. In our scenario, the 
coefficient of this 2PN correction is 0(1) because I is re- 
lated to the total mass M as l/GM ~ 1. This correction 
affects the relation between the energy and the orbital 
angular velocity of the binary and hence the frequency 
evolution of the binary is modified. As for the energy loss 
rate from the binary, one can show that the correction is 
of higher order. 1 

Fourier components of gravitational waves from an in- 
spiralling binary in a quasi-circular orbit are given by 

M/) = |(GM c ) 5 / 6 /- 7/6 exp(i*(/)), 

where Q is the factor depending on the direction and the 
orientation of the binary and the detector antenna pat- 
tern. Keeping the second order term in e, the phase ^f(f) 
can be written as 



*(/) = 2nft c 
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where 77 = (J./M, t c is the coalescence time, (3 denotes 
the spin parameter, and <p c determines the phase. The 
KK mode correction appears in the last term in j(f). 
Here we have neglected the spin dependent term at the 
2PN order. One may think that the KK modes contri- 
bution is indistinguishable from this effect. However, for 
binaries whose spin parameter at the 1.5PN order (3 is 
measured to be small, one can neglect the spin effect at 
the 2PN order. Hence, the effect of extradimensions can 
be observed as an anomalous "2PN" correction. 
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FIG. 1. Sensitivity curves h„(f) = Sl' 2 {f) for LIGO-II 
and EURO with different specifications, (A,B,C) fk = 10 3 Hz, 
fk — 6.5 x 10 3 Hz, and the xylophone type, respectively. 
Disks and circles denote the wave strength h s for the near- 
est PBH SULCO binary with mass M B h = KT 7 M for the 
observational time AT = 10, lyr, respectively. Dotted lines 
and dashed lines represent the evolution of the amplitude 
Q(GM c ) 5/6 f- 7/6 /D min for the nearest PBH SULCO binary 
with Mbh = 1CT 7 M , Mbh = lCr 8 M , respectively. The 
endpoint of each line in the left-upper side corresponds to the 
remaining time T = lOyr before coalescence. 

Let us evaluate how accurately we can determine the 
AdS radius I, using third generation detectors such as 
EURO. We set f k = 6.5 x 10 3 Hz in the following analy- 
sis. We assume that a PBH SULCO binary with lunar 
mass Mi = M 2 = 10~ 7 M Q is in a quasi-circular orbit. 
With AdS radius I = M~ 2 M = 0.3mm, a binary that 
emits gravitational waves at a frequency / = 2.5 x 10 3 Hz 
coalesces in 10 yr. The distance to the nearest binary is 
estimated as D m in ~ 3 x 10 2 pc, which yields the ampli- 
tude h = Q(GMc) 5/6 D-\f- 2 ^ - 1 x 10- 24 . 

In the matched filter analysis up to 2PN order for 
the phase ^(/), we marginalize seven parameters {Ai} = 
{In A, t c ,<p c , lnM c , In fi, (3, In I}. The error of each param- 
eter can be estimated by a quadratic fitting of the likeli- 
hood function assuming uncorrelated Gaussian noise [15]. 
Using the EURO detector, we find that the nearest PBH 
SULCO binary is expected to be detected with a signal- 
to-noise ratio S/N = 3. The relative error in I is 
Al/l = 3x 10" 1 whereas At c = 2 x 10~ 4 sec, A(j) c = 
5 x 10 2 rad, AM C /M C = 7 x 10" 11 , Afi/fi = 7x 10~ 5 and 
A/3 = 4 x 10~ 2 . Although the S/N is not sufficiently 
large, the AdS radius I can be determined within <~ 30 



1 According to the AdS/CFT correspondence, the energy loss owing to geometrical particle creation in the RSII model is 
30-rrl 2 /G times as large as that for a single conformally coupled scalar field [13]. Using the result for a scalar field [14], one can 
easily find that the correction is suppressed by a factor f 2 £ 2 , and hence is of 3PN order. 
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percent error. To gain a large S/N value, we need a de- 
tector that is highly sensitive at high frequency bands 
(>lkHz), say, a "xylophone" type detector [11] that con- 
sists of several narrow-band interferometers. 
Summary and Discussions) In this letter, we discussed 
observability of gravitational waves from the nearest co- 
alescing PBH SULCO binary. If PBH SULCOs with 
mass Mbh ~ 10~ 7 M Q constitute the cold dark matter, 
the third generation interferometer, EURO, has sufficient 
sensitivity to confirm their existence. For the case that 
the PBH formation is related to the braneworld scenario, 
we further discussed observability of the imprint of ex- 
tradimensions. We found that the sensitivity of EURO is 
marginal to discriminate the RSII model. However, once 
we discover the PBH SULCO binaries, we would be able 
to construct a new detector that is suitable for probing 
the extradimensions. 

Although we have investigated the RSII model as an 
example, our estimate of observability is not very restric- 
tive to this model. In the RSII model, the deviation from 
general relativity starts at the 2PN order. In other mod- 
els the deviation may appear at lower PN order. Even 
in that case, the correlation with other parameters such 
as M c , ji and (3 renders new physics unveiled up to the 
2PN order. Therefore, it is plausible to conclude that our 
result is generic to a rather wide class of models. 

Because the number of PBH SULCOs in the dark halo 
is huge, one might expect that the radiation from gas 
accretion to the PBHs is detectable [16]. However, as 
we shall see, the chance of detection is hopelessly small. 
Let us consider that a lunar mass Mbh = 1O _7 M0 
BH is moving with velocity V with respect to the sur- 
rounding gas. The accretion rate from the surround- 
ing gas can be approximated as m = Mbh/ MEdd ~ 
5 x 10- 4 (M Bff /M Q )(n 3 /10 2 cm- 3 )(10kms-7U) 3 where 
n g is the number density of the surrounding gas, and 
the MEdd is the Eddington accretion rate with unit ef- 
ficiency. Adopting the thin-disk model [17], the to- 
tal luminosity from the disk is estimated, at most as 
L <~ GMbhMbh I '(2J? TO m) where R m in is the innermost 
accretion radius. Assuming black-body radiation from 
the disk, m = 1 and R rnin = 3R Sc hw where Rschw 
is the Schwarzschild radius, the temperature becomes 
T~lx IO^Mbh/Mq)- 1 / 4 . Thus smaller mass BHs 
have a higher temperature. Now consider lunar mass 
PBHs moving in the Orion Nebulae at 400 pc away. As- 
suming that the dark matter density in the solar neigh- 
borhood is O.OO79M pc~ 3 [10], the total number of lunar 
mass BHs is just <~ 8 x 10 7 in the molecular cloud if its ex- 
tention is 10 3 pc 3 . Provided that the velocity distribution 
of the PBHs is Maxwellian, and the number density of the 
molecular gas is n g = 10 2 cm~ 3 , we find that the expected 



total luminosity is L^2 x 10 24 ergs/s. and the peak fre- 
quency is / ~ 10 18 Hz. On the other hand, the sensitivity 
of the X-ray telescopes such as the Advanced X-Ray As- 
trophysics Facility (AXAF) is about ~ 10 29 ergs/s for a 
source at 400pc at / <~ 10 18 Hz. Thus we cannot expect 
any X-ray visible BHs in the Orion molecular cloud. 

Based on the AdS/CFT correspondence, some authors 
have conjectured that the evaporation rate of BHs in 
the RSII model is greatly increased [18,19]. The life- 
time of massive BHs in the RSII model is estimated as 
r - 10 5 (M BH /10M Q ) 3 (0.1mm/Z) 2 yr for a BH with a 

mass Mbh 3> IM^ [18,19]. Thus, most astrophysical 
BHs with moderate mass 10M Q < M B h < 10 3 M Q would 
have evaporated by now, leaving remnants with a small 
mass ~ IMp. If this should be the case, we would be able 
to observe gravitational waves from BH SULCO binaries 
that are not primordial ones, though the observational 
evidence of stability of stellar mass BHs would give a 
stringent constraint on the size of the AdS radius I. 
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